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Abstract. A theory for the electromagnetic response of type-II superconductors close 
beyond the critical state is presented. Our formulation relies on general physical 
principles applied to the superconductor as a thermodynamic system. Metastable 
equilibrium critical states, externally driven steady solutions, and transient relaxation 
are altogether described in terms of free energy and entropy production. This approach 
allows a consistent macroscopic statement that incorporates the intricate vortex 
dynamic effects, revealed in non-idealized experimental configurations. Magnetically 
anisotropic critical currents and flux stirring resistivities are straightforwardly included 
in three dimensional scenarios. 

Starting from a variational form of our postulate, a numerical implementation 
for practical configurations is shown. In particular, several results are provided for 
the infinite strip geometry: voltage generation in multicomponent experiments, and 
magnetic relaxation towards the critical state under applied field and transport current. 
Explicitly, we show that for a given set of external conditions, the well established 
critical states may be utterly obtained as diffusive final profiles. 
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1. Introduction 

Over a half century, the critical state model [1] (GSM) has been an essential 
interpretative tool for the investigation of the macroscopic magnetic properties of type-II 
superconductors. This accomplishment probably relies on a clear physical background 
together with an ostensible mathematical simplicity, that have enabled a huge number 
of utilizations both in material characterization as well as in more fundamental studies. 

In brief, the GSM postulates that the biased magnetic response of a type-II 
superconductor is defined in terms of a series of metastable equilibrium configurations, 
each characterized by a current density distribution J(r). In the simplest cases, J(r) may 
be obtained by integration of Ampere's law given by dH/dx = ±Jc , 0. Jc, the so-called 
critical current density is the single material parameter of the theory, and characterizes 
the balance equation between magnetic and intrinsic pinning forces: J x B = Fp. The 
transition between different configurations is assumed to take place instantaneously, 
which means that external field variations occur slowly enough as compared to the scale 
established by the material response (Text ^ Tmat)- Thus, although it is well known that 
magnetic flux penetration in the presence of pinning forces happens as an avalanching 
process [2] when the threshold condition is exceeded (i.e.: J > Jc), one argues that 
magnetic diffusion is so fast that the superconductor settles a negligible time in the 
intermediate resistive states. 

In many instances, the GSM approximation is justified, but especially for the case of 
high Tc superconductors one can meet practical situations for which relaxation transients 
towards equilibrium configurations have to be considered. A very relevant case is the 
fault current limiter that, by construction operates driving the superconductor beyond 
the resistive transition [3l|4]. The precise knowledge of the magnetic diffusion processes 
that occur is obviously necessary for the design of such devices. In particular, the 
essential evaluation of energy losses may be seriously tampered if one neglects the effect 
of the sample becoming resistive. On the other hand, and more to the side of basic 
physics, one has to recall that, experimentally, the value of the critical parameter Jc 
is often obtained from transport measurements that are based on some threshold value 
for the voltage detected when the sample goes resistive (typically IfiV/cm). In a recent 
work about the origin of the dissipation mechanisms that operate for superconductors 
in the vortex state, it has been argued that, although designed for obtaining the critical 
parameters, resistive measurements could be ineluctably recording properties of the 
resistive state [5]. This circumstance is relevant for HTS samples that display more 
gradual current voltage transitions as compared to the sharp behavior for conventional 
superconductors. Finally, it should also be mentioned that from the fundamental point 
of view, embedding the GSM approximation in a time dependent theory that allows to 
derive it as a limiting case is by itself a desirable objective. 

In this work we put forward a minimal formulation that upgrades the critical state 
theory so as to enable the inclusion of overcritical behavior. The basic picture of our 
statement is as follows. Subject to an external action, the state of the superconductor 
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will drift until new steady conditions are met. If allowed, the steady state will be an 
equilibrium phase, characterized by a conventional force balance equation, with critical 
behavior for the current density. On the other hand, dissipative steady states and 
excursions from one critical state to the other will be described under a first order 
(linear) approximation for the underlying driving forces. Following the spirit of Bean's 
model, although describing phenomena whose ultimate nature would at least require 
a mesoscopic scale, we pursue a theory based on macroscopic variables. Thus, our 
basic assumption is that the superconductor behaves as a thermodynamic system, with 
the coarse-grained current density as the state variable, and postulate that dynamic 
properties may be predicted on the basis of three factors: 

(i) magnetic forces 
(i) material pinning forces 
(iii) irreversible thermal forces. 

The latter group will allow to introduce the physical effects related to dissipation 
induced by overcritical current flow. In order to avoid microscopical modeling, a power 
series expansion argumentation will be used for introducing such effects. Outstandingly, 
recalling universal properties as entropy increase and single valuedness of the physical 
observables, one can fairly determine the kind of models to be used, even in complex 
scenarios such as flux cutting environments or 3D modeling. 

It must be clarified that, being interested in the description of small excursions 
beyond the critical state, long term effects, such as thermal relaxation to true equilibrium 
are not considered here. Thus, our theory describes the behavior of superconducting 
material with a characteristic a time very short as compared to thermal creep relaxation, 
but non negligible against external source variations, i.e., Tmat ^ ^cxt ^ Tthermai- 
Therefore, magnetic diffusion times are small but not negligible against external sources 
variations, while thermal creep (of the order of hours or days) is ignored, determining 
a well defined metastable equilibrium. Furthermore, isothermal conditions will be 
assumed and, thus, any energy transfer from the electromagnetic sector would be 
immediately absorbed by the surrounding thermal bath. 

The article is organized as follows. Section [2] is devoted to deliver the formal 
details of our theory. First (Sec j2.ip we recall some mechanical concepts (fields, forces, 
Drude's model) that allow a rudimentary approach to the problem of dissipation in 
normal conducting systems. Then, (Sec J2.2l) some basic thermodynamic background is 
introduced, with the aim of paving the way for the generalization to superconductors. 
As a central result of our work, in 12.2.21 we introduce the idea of dissipation function 
and entropy generation for type-II superconductors driven out of equilibrium. We will 
show that, concomitant with a complex structure for the metastable equilibrium states, 
an acceptable theory of resistive losses has to fulfill some consistency requirements, 
that will be used to put restrictions on the possible material laws to be used (E(J) 
in particular). In order to ease the practical implementation, a variational form of 
our theory is issued. The second part of our paper (SecE]) illustrates the application 
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of the above concepts to practical situations. Mainly focused on the infinite strip 
geometry (quasi-lD configuration), we show that our basic equation allows to obtain 
the electromagnetic quantities either in equilibrium, steady states or during transient 
processes. It is explicitly shown that, for a given set of external conditions, relaxation 
eventually leads to the well known corresponding critical states. Excellent comparison 
with analytical results, when available, is displayed. The final section (SecH]) summarizes 
our results and contains a brief discussion of possible applications and extensions of our 
work in the area of type-II superconductivity. 

2. Resistive losses in hard superconductors 



2.1. Coarse grained modeling: fields, forces and E{J) law 

Let us start by recalling some details about the classical description of electrical 
conduction in normal metals, viz., the Drude model. As we will see, this simple scheme 
may be illuminating for the issue of a minimal model of the resistive behavior of hard 
superconductors. 

2.1.1. Normal metals 

The simplest description of electric current in normal conductors (Drude model) 
is built from the classical dynamical equation of the charge carriers subject to both 
an external electric field and a phenomenological drag force Fdrag standing for the 
interaction with the molecular environment and other charges, i.e.: 



Customarily, Fdrag is taken to first order, i.e.: Fdrag ~ — ('^eVe/rtr) and it is a 
characteristic of the material through the time constant r^r- It is important to recall that 
all the quantities in the above formulae have to be interpreted as mesoscopic averages 
both in spatial regions and time intervals. This entails smoothing (by randomizing) the 
macroscopically irrelevant fluctuations of the microscopic level. In a thermodynamic 
language, one speaks about a thermalized electron-lattice system. In spite of its 
simplicity, the model catches the basic behavior observed in an overwhelming amount 
of experiments. 

It is of particular interest to observe that, unless for very high excitation frequencies 
(in the range of 1 THz for typical metals), further simplification is possible. First, a 
quasisteady state approximation may be used, that implies to neglect the time derivative 
(electron inertia) to the left hand side of Eq.(IT]). Then, the non-dispersive form of Ohm's 
law arises, i.e.: J = (?T,e^rtr/?TieE) = ctqE with n the number of conduction electrons 
per unit volume. Also, related to the use of good conductors, as it is customary in 
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conventional electric machines, charge accumulation within the sample may be neglected 
(i.e.: dpq/dt ^ 0) and one may use Ampere's law in its magnetoquasistatic form 
V ■ J = 0. In practice, when studying the evolution induced by some process of the 
external excitation, this means that one starts by calculating the dominant magnetic 
fields from the diffusion equation, that for latter convenience will be expressed in terms 
of the material resistivity po = V'^o 



Eventually, E is obtained from Faraday's law. Notice that stationary solutions of this 
equation verify the condition V x J = 0, that is to say: persistent current loops within 
the sample are excluded for a normal metal. 

2.1.2. Hard superconductors 

Microscopically, the nature of resistive losses in flux penetrated type-II 
superconductors is quite more complex than the scattering of normal electrons, which 
is behind the above summarized Drude's model. Nevertheless, one can take benefit 
of the conceptual aspects introduced, that may be straightforwardly applied at the 
phenomenological level. Thus, it is well known that under the action of a transport 
current and a perpendicular magnetic field, these materials experience a resistive 
transition, also characterized by a linear E( J) relation, now in terms of the so-called flux- 
fiow resistivity pfj. This behavior is well understood at the mesoscopic level. Abrikosov 
vortices are driven by the Lorentz-like force J x B and their normal cores contribute to 
the transport current as a normal channel. Then, in the absence of additional effects, 
one could formulate the behavior of the superconductor in the same terms described 
above for a normal metal, just replacing po by pg. 

However, the detrimental flux flow behavior has been customarily attacked by 
a number of pinning strategies. Basically, they rely on the idea of introducing a 
restoring force on the vortices, so as to keep them in equilibrium positions while a 
transport current flows along the system (J x B = Fp). As pinning forces are bounded 
{Fp < Fp^rnax), there will be a threshold value for the current density that can flow 
lossless, the so-called critical current density Jc- As it was said before, in this framework, 
the calculations leading to the evaluation of magnetostatic equilibrium properties may 
be rather simple. Thus, Bean's postulate states that flux penetrated regions will be 
characterized by the threshold (critical) conditions, i.e.: \J\ = Jc- Following this, in ID 
problems one just has to integrate Ampere's law in the form dH/dx = ±Jc, supplied 
by specific boundary conditions for the magnetic field. Nevertheless, the investigation 
of the transient processes involving the appearance of resistivity may be quite complex, 
even for the simplest geometries. The reason is that now, effective drag forces only occur 
for currents circulating with a density beyond the critical value Jc- In other words, one 




(2) 
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Figure 1. (Color online) Piecewise linear E{J) law for a type-II superconductor. The 
material parameter Jc defines the boundary Tc ^ { — Jc,Jc} of the lossless behavior 
(region A = [—Jc,Jc\ on the horizontal axis and detailed as an inset). Ohm's law 
corresponds to Jc = => A = 0. The arrows indicate the effect of introducing pinning 
forces, and thus progressively increasing Jc- 



has to deal with a material law of the kind E = p{J)J where 
, / if |J| < J, 

As a consequence of the non-linearity introduced by the piece-wise constant 
behavior of p(J), one has to deal with noticeable difficulties, as for instance, the form 
taken by the diffusion equation (compare to Eq.([2])) 

(^/.o^ - p(J)V2^ H = (V X H) X Vp(J) (4) 

Difficulties for solving this equation by analytical methods have been reported in 
Refs.( OH]) that basically suggest to apply ad hoc approximations when justified by 
the set of experimental data under consideration. On the other hand, we stress that, in 
fact, the limiting Bean's approximation means to avoid such equation under the following 
thoughts. If the actual material law in Eq. ([3]) is such that one can speak about a sharp 
transition, i.e.: takes elevated values as compared to the experimental parameter 
E / J , one can approximate pfj — oo. Then, the magnetic diffusion time (xp oc 1/p) may 
be neglected and externally induced evolutions may be accurately described as a series 
of equilibrium critical states. 

Fig{T] displays the main features of the above mentioned modeling with a finite 
value of pff for the case of ID systems. Let us analyze the physical significance of such 
a material law. 
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(i) Note first that superconducting lossless transport is allowed by the condition E = 
within the region A = [— Jc, Jc]- The superconducting range A increases with the 
value of Jc, starting from in very low pinning materials, that behave like "normal 
conductors" (see arrows in Fig{T]). 

(ii) Persistent shielding is also possible, subsequent to electromagnetic induction 
because E can go to zero and J take a local structure that allows V x J 7^ 
when diffusion stops. In particular, this includes the critical state solutions with J 
taking the values ±Jc or within the sample, and forming screening loops. 

(iii) One can visualize Bean's model as a limiting case of this material law when there 
is an arbitrarily high slope of E{J) for J > Jc- 

(iv) The conceptual implications of FigiUmay be straightforwardly generalized to higher 
dimensions. In particular, A may be a region of the J-space and Tc its boundary. 
Thus, for the case of high-T^ superconductors it has been recently shown [5] that 
a meaningful selection would be to take A as an elliptic region, within the plane 
defined by the components of J parallel and perpendicular to the local magnetic 
field. Other possible regions, related to specific materials, will be commented on. 

(v) The material law for a type-II superconductor depends not only on the intrinsic 
parameter ps, but is strongly affected by extrinsic quantities. In particular, the 
ambient magnetic field will play a prominent role for the definition of the region A. 
Thus, already in the simple ID cases, the dependence Jc{H) is usually an important 
concern. On the other hand, as mentioned above, the relative orientation of the 
vectors J and H is crucial for the material law in higher dimensions. 

Performing calculations within the above described material law and its natural 
extensions to 3D is by no means a simple task. Below, we suggest a first approximation 
to the problem that relies on thermodynamic concepts applied either close to equilibrium 
or close to steady states. This will provide a useful framework for developing numerical 
tools that allow to analyze the resistive transition in applied superconductivity. 

2.2. Thermodynamic background 

Non-equilibrium thermodynamics is still an area under development. One of the main 
fields of interest, i.e.: the study of steady states will guide us through the analysis of 
electrical conduction. In fact, being concentrated on quasisteady systems, the extension 
to irreversible transient processes (magnetic diffusion) will be straightforward. As an 
example of what can be done, we recall [7j that based on the law of increasing entropy 
and also on Onsager reciprocal relations, one can show that in normal conductors, 
the conductivity tensor must be positive definite and symmetric. This generalizes 
ID Ohm's law. Obviously, such properties could never be deduced from the pure 
electromechanical principles used in the previous paragraph. Thus, having the aim of 
a theory that generalizes Sec J2.1.2l for type-II superconductors, in the forthcoming, we 
present a minimal conceptual scenario that allows to host the main physical properties 
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of these materials. As before, we will start introducing the main ideas by taking benefit 
of our knowledge of the simpler normal conductors. 

2.2.1. Entropy and dissipation function in normal metals 

Let us start by recalling some definitions. Being interested in local properties of 
the conductor, we consider the mesoscopic entropy function S as an average that has 
smoothed the statistical microscopic fluctuations. In fact, the spatial dependence comes 
through the current density, that will be our state variable: iS[J(r)]. Here, r denotes 
the position of a region of mesoscopic size, and nonlocal correlations are neglected. In 
the absence of an external action the system settles at the entropy maximum, i.e.: the 
stable equilibrium point J = 0. Physically, an overwhelming amount of all possible 
microstates corresponds to macrostates with J ~ 0, and then 5(J) < 5(J = 0) = 5eq. 
Statistical fluctuations around such point become negligible with increasing size of the 
subsystem. 

Out of equilibrium, but not far from it, the same statistical mechanism that 
suppresses fluctuations will operate. Thus, after a displacement caused by an external 
agent (electric field in our case), a thermodynamic "restoring force" Fdrag in Sec J2.1.11 
drives the system along increasing entropy until 5eq is reached again. F^rag is responsible 
for the energy losses that may be expressed in terms of the so-called dissipation function 
[Sj J-" = (l/2)Fdrag ■ Ve, that is a measure of the amount of heat generated per unit time. 
Macroscopically, the amount of heat generated is given by [9] 

dQ = TdS = 2Tdt . (5) 

As stated before, isothermal conditions will be assumed, and we will use the relation 
S = 2J^/T. Now, thermodynamics enters through the so-called "principle of minimum 
entropy production" [lOlE] that may be expressed as 

The steady state of a system is that state in which the rate of entropy production 
has the minimum value consistent with the external constraints. 

Notice that, based on this, the equilibrium state takes its natural place when there 
are no constraints on the system, because as a consequence of the second law, it reaches 
the maximum entropy, and thus the absolute minimum of entropy production, i.e.: zero. 

In order to see how the above principle operates, let us assume that some external 
action (electric field) drives the system out of equilibrium. If the constraint is removed, 
a transient towards equilibrium (Scq) will occur. Around this point, the function will 
allow a quadratic expansion of the kind 

^^ij^J.l^^^J,. (6) 

with fi*-' the components of a symmetric positive definite tensor. It is apparent that, 
in the absence of constraints, minimizing J-" leads to the condition J = 0, and that 
the thermodjTiamic field driving the system towards equilibrium may be obtained 
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as Ether = — VjJ-" = —Q3. Its lines of force (maximal slope of J-") in the J-space 
correspond to the most probable macroscopic evolution (faster increase of entropy). In 
steady situations, this field balances the applied electric field that verifies Ohm's law 

(-Ether = E = fij). 

On the other hand, steady states out of equilibrium will be characterized by a 
constrained minimization problem, in which J-' has to be augmented by some Lagrange 
multiplier. For instance, as indicated by Landau [7], the steady current distribution 
within a normal conductor may be obtained by minimizing the volume integral of 
J-" + A VJ as corresponds to the stationarity condition V J = 0. 

As a main conclusion from the above paragraphs, we resolve that the behavior of a 
normal conductor close to equilibrium may be formulated in terms of a quadratic tensor 
(resistivity) that is positive definite and symmetric as dictated by thermodynamics. 
This obviously generalizes the one dimensional results in Sec J2.1.11 An additional 
observation is that just recalling the definition of the gradient function one can show that 
the electric field lines will be perpendicular to the constant level sets of the dissipation 
function J-". Finally, variational principles allowing to afford complicated calculations 
as those related to the diffusion towards the equilibrium point (J = 0) or characterizing 
steady sates may be stated in this realm. 

Below, we show that the ideas in this section may be exported to the case of type-II 
superconductors. However, already related to the generalization of the one dimensional 
E{J) law in Fig{Tl some technical difficulties arise, and further analysis is required. 

2.2.2. Entropy and dissipation function in hard superconductors 

In the light of Sec J2.1.2l the concept of thermodynamic equilibrium has to be 
reconsidered for pinned superconductors. Thus, even for the simplest one dimensional 
configurations, one finds a vanishing drag force not just for a point (J = 0) but within a 
whole segment (region A) that becomes a surface, or even a volume in the J— space for 
more complex scenarios. Then, the entropy of the system will be maximal within a full 
set of values Jc € A and the determination of the state from thermodynamic arguments 
seems flawed by ambiguity. In fact, a complex multiplicity of possible states would 
seem admissible. Against this, one could argue that, in fact, either static or stationary 
configurations always occur subsequent to some diffusion process accompanied by drag 
forces (electric fields) and that "integration" along some path within the J-space that 
connects the overcritical region and the eventual critical point Jc can be performed. 
However, the following mathematical quiz arises: how does one expand a function J-" 
towards the resistive behavior when a "starting region", instead of a starting point, is 
given? 

Once more, a very general argument, uniqueness, will help. Let us consider the 
kind of elliptical region mentioned before, that is physically meaningful as related to 
dissipation mechanisms in either Jy or J± relative to the local magnetic field orientation. 
As depicted in Figj2l irreversible energy dissipation corresponding to some overcritical 
point J ^ A may be uniquely quantified by means of d, i.e.: the minimum distance 
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from the point J to the boundary Tc- Geometrically, this means that expansion is done 
perpendicular to Tc- Analytically, one has to find some expression that allows to obtain 
the critical point Jc for each value of J. In the case of the elliptic region considered 
here, this entails to solve a quartic equation. Of special mention is that meaningful 
cases [12] as the isotropic model (A is a circle) and the Double Critical State Model (A 
is a rectangle) produce trivial conditions for the determination of the point Jc- 




Figure 2. (Color online) Dissipation function for a type-II superconductor close to the 
critical state. The horizontal plane is defined by the current density vector components 
either along (J||) or perpendicular {J±) to the local magnetic field. For a given value of 
the current density vector ( J|| , Jj_), is defined through the distance d to the critical 
region boundary Tc- The related critical point Jc = (J|j, J^) is also shown. 

From the physical point of view, the above mathematical conditions imply that 
for small perturbations around equilibrium, the induced electric fields will keep 
perpendicular to the levels of constant dissipation J-". This extends our previous result 
[12] that just "on surface" (critical state) E is normal to the boundary of the region A. At 
least for small perturbations, a consistent theory that assumes the existence of a critical 
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region with boundary may be obtained by using E ± and prolong this condition 
towards the resistive state. We call the readers' attention that this general property is 
not restricted to the case of a region A defined in terms of J\\ and J±. This would just 
fit the case of a homogeneous type-II superconductor, whose anisotropy only depends 
on the local magnetic field. Other physical scenarios, such as microstructure induced 
anisotropy are apparently tractable, just by using the required axes. For instance, one 
could use J"'' and J'^ in the case of uniaxial symmetry for the current fiow. A would be a 
cylinder with its axis parallel to the c direction. Also, one could introduce a dependence 
of the critical current density on the orientation of the local magnetic filed relative to 
some axis in the material along which pinning is favored, i.e.: define Fc by a function 
Jc{<y) = Jc(H. ■ n) = Jc(cosa). In fact, the critical state solutions of such statements 
were already investigated in Ref. [13] and the extension to the fiux fiow regime may be 
done in the above terms. 

Before pursuing a specific mathematical statement of the above concepts, some 
clarifying words on the physical system are still due. On the one side, we have 
invoked the molecular thermal forces so as to describe the dissipative behavior beyond 
equilibrium. On the other side, criticality relates to a multiple number of possible 
"equilibrium" states for the system. The full framework is as follows. Indeed, thermal 
forces acting on the long term would eventually drive the system to a true (and unique) 
thermodynamic equilibrium state with J = 0. However, as we are focused on small 
perturbations around the critical state, we just consider the fast relaxation towards 
a given metastable configuration determined by the applied fields and pinning forces. 
This corresponds to recalling that, in hard superconductors, much creep occurs in the 
first instances when the fiux density gradient (basically, our Jc parameter) is high |14j . 
Subsequent flux homogenization by flux hopping can be neglected in the magnetization 
processes under consideration as far as its typical time would be much beyond the 
experimental period. 

2.2.3. Thermodynamically admissible EfJj laws: mathematical issues 

The dissipation function F must be deflned, positive and obviously single valued, 
outside the region A of (sub) critical current densities, in the J vector space 

J" : IRVA ^ IR+ J^(J) > (7) 

A Taylor expansion around the boundary of critical currents F^ will encode the main 
properties of the electromagnetic behavior close beyond the critical state, i.e., the E(J) 
constitutive law. 

To start with, we will consider spatially isotropic samples. Thus, the surface Fc will 
be axially symmetric around the local H direction, as well as symmetric under reflection 
through the plane normal to H. A is also supposed convex. With this symmetry, the 
local electric fleld takes value in the plane deflned by H and J. Hereafter, we restrict the 
analysis to this plane, and use cartesian coordinates (^1,^2), (-^1,-^2) and (ifi,0). We 
note that the cartesian coordinates are nothing but the physical directions (||, J-) deflned 
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Figure 3. (Color online) Schematics of the current density vectors in the over-critical 
state problem (as in Figl2]). J stands for the actual current density, Jc for the nearest 
critical current density vector, and AJ for their difference. Related angles are also 
defined. 



in Figs|2] and [31 However, in order to avoid awkward notation, within this section we 
use (1,2) = (II, ±). 

Fc is then a closed curve surrounding the convex region of subcritical currents. 
In polar coordinates (Ji, J2) = (Jcos(^), Jsin(^)), and is determined by a function 

Jc{e). 

Let us see how the concept of perpendicular expansion introduced above arises. 
In principle, there is no particular point at Fc from which to perform the Taylor 
expansion. Let us then consider an arbitrary ji^^ G Fc, and denote by T(J; ji^^) the 
Taylor expansion of J-" ^ T(J; ji^^) around this point Jc^\ It is reasonable to take a 
point ji^^ close to the value J of interest, in particular, nearest point to J on the critical 
surface. However, a well behaved function J-' should admit compatible approximations 
from nearby points, say, J-" ~ T(J; Jc^^). Therefore 

T(Jo;J«)^T(Jo;jf)^T(Jo;JeW) deT = (8) 

i.e., when fixing a value Jq, the derivative of T(Jo; Jc) along the curve F^ must vanish. 
Now, for our purposes it suffices to include second and third order terms 

T(J; J«) = l[pn( Ji - Jil^r + P22{J2 - 4?)^ + 
+ 2p^2iJl-Jil^)iJ2-J^J)] + 

+ J[7iii(^i - + 37ii2(Ji - Ji'JYiJ2 - + 

+ 37122(^1 - J[1^){J2 - 4J? + 7222(^2 - 
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with 

From the condition dgT = we get, to first order 

(pn,Pi2)-T(^) = (pi2,P22)-T(e) = (10) 

so that T = dg{Jc{9) cos{9), Jc{0) sin(6')), the tangent vector to Fc, belongs to the kernel 
of the resistivity tensor, which in physical terms means that no resistive losses occur for 
paths along F^ 

n=( ] (11) 

y Pl2 P22 J 

On the other hand, to second order 



(7111,7112) • T{9) = - dgpu 
(7112,7122) ■ T{9) = - depi2 
(7122,7222) ■ T(6') = - dep22 

Note that for the strict second order Taylor polynomial (74^^ = 0), the compatibility 
condition implies that Vt is constant, but this is no longer true when third order terms 
are considered. 

Now, if we use the vector basis {t,n}, i.e.: unit tangent and normal vectors at a 
given point of Fc, the tensor VL becomes 



n=K (12) 




Thus, the Hessian of J-" is positive semidefinite, with t null eigenvector, and pn the 
eigenvalue of fi. 

Remarkably, one can maintain the simplicity of the second order approach while 
allowing anisotropic ^2(6'), by performing the Taylor expansion as follows; from each 
point Jc G Fc we compute the Taylor expansion "exclusively" along the normal halfline 
Jc + AJn. In such a way, and taking into account the convexity of the base curve, there 
is a univocal correspondence of each J with a particular Taylor expansion, and the 
whole region IR^/Ac is covered. The second order approach to the dissipation function 
becomes 

= lp(a)(AJ)2 J = Jc(a) + AJn(a) (13) 

The correspondence J = (J, 6) ^ (A J, a) is the geometric condition of nearest point 
of the curve Fc to J, and it obviously depends on the specific Jc{oi) function. Note the 
distinction between 9 and a, the angular coordinates of J and Jc (Fig|3]). 
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Eventually, the electric field (E = Vj J-") is given by the polar coordinate expression 
E = ajJ^J + -deJ^e = E{3 + eJ . (14) 

For practical purposes, an intrinsic coordinate system representation will be useful. 
Thus, a better adapted expression in polar-like coordinates {AJ, 0}, with A J the 
distance of J to Fc, and the angle between AJ and the H axis reads 

E = aAj^AJ + i9^^0. (15) 
Itj) 

Here is the length of d^J = d^Jc + A JS^n, and AJ = n. Then, one has 

E = p(0)(AJ)n + -La^p(AJ)20 (16) 

Notice that, in case of a constant p, the electric field is always parallel to fi. However, 
as it will be analyzed below, in connection with voltage-current experiments, the ratio 
E\\ and E± changes with the separation of the working point and the critical region 
A J = J — Jc (see section [3^ . 




Figure 4. Schematics of a simple irreversible system: a charge carrier subjected to 
the action of an electric field and a simultaneous viscous force. 



2.2.4- Variational formulation of the conduction problem 

It has been argued that classical mechanics and its methods do not provide a 
complete framework for the analysis of electrical conduction. However, having exploited 
the consequences of the second law of thermodynamics, one can reconsider the problem. 
In particular, below we show that by including the concept of admissible dissipation 
function J-" one can issue a variational formulation for the magnetic diffusion problem 
between steady conduction states. In this unified description that encompasses 

normal conductors and type-II superconductors is presented. 

Just with illuminating purpose we start by considering a ID problem with a charged 
particle subject to an electric field E with associated potential energy U. Recall that 
the Lagrangian formulation of Hamilton's principle is as follows 



S = / L(x, v) dt 



U{x) 



dt 



_ , ddL dL dv dU ^ 

MmS ^ __-_ = ^ m— = -— = E 

dt ov ox dt ox 



cons 
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where Fcons stands for the conservative force Fcons = qE. 

Consider now that a viscous drag acts on the particle (as depicted in Fig |4] for the 
case of an electron). Notice that a minimum principle leading to the sound equations 
of motion can still be formulated for a modified Lagrangian as shown below 

S = j L{x,v,t)dt with L = L + I't = L + ^hv^t 

r ddl dl ^ dv dU ^ ^ 

MmS ^ - =0 ^ m— ~ hv = F^ons + F^ 

dt ov ox dt ox 

What has been done in deriving such formulation is to neglect variations of 
the viscous force F^ within the interval of time considered. Then, the suggested 
approximation will be valid if minimization is applied "iteratively" with intervals of 
duration much less than the characteristic time t = h/m. This relates to the so-called 
Adiabatic Hypothesis used in other physical disciplines: 

// energy, though not conserved, varies slowly according to some parameter, then 
one is allowed to assume a kind of isolated system within small enough intervals of 
the temporal evolution. 

An eventual calculation will help us to gain some more insight. Let us use the classical 
mechanics expression for calculating the energy of the particle in terms of the modified 
Lagrangian. Upon neglecting the term AF^Ax the energy within a given interval is 

dL f mv^ 1 ^ * /X 

that is to say, we are effectively dealing with an isolated system which stores an average 
energy of one half of the full loss F^vAt. 

Let us now see how the above arguments may be exported to the problem of 
electrical conduction. Within the quasisteady approximation, the incremental time- 
averaged field Lagrangian for conducting materials reads 

(L) = ^[ ||Hn+i-H„||W+ / AtJ^dV={[ CdV) (18) 
2 jr3 Jvoi Jr^ 

with J-" the dissipation function introduced in Sec J2.2[ Here, Hn means the magnetic 
field distribution at a given time instant given by tn and At = tn+i — tn- 

For the case of a normal conductor, it is relatively simple to show that the Euler- 
Lagrange equations lead to the desired diffusion equation (Eq.([2])). Thus, if one assumes 
a diagonal resistivity matrix, i.e.: = po^*"' and replaces J-" by its value, it follows 



dC 



E 



d 



dC 



dxj d ( 9Hn+ 1 / dxj ) 



Hn+l — Hn 

At 



PoV X V X Hn+l = 



(19) 
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i.e.: the time discretized form of Eq.(|2]). 

Outstandingly, the interest of Eq. ffTSj) is that it gives way to the possibihty of 
applying direct numerical minimization methods. In particular, this allows to deal with 
non-simple forms of J-" as required for the investigation of type-II superconductors (recall 



From the mathematical point of view, one can distinguish two kinds of problems 
as related to the minimization of the functional in Eq.f llSp : (a) strictly variational 
structures, when either the first or the second term may be neglected, and (b) a quasi- 
variational structure when both are relevant [15]. Physically, one would speak about: 

(i) Equilibrium-like processes, leading to the critical state (the dissipation function may 
be neglected and one basically obtains the configuration that balances magnetic and 
pinning forces). 

(ii) Steady states within a dissipative regime (the magnetic inertia AH may be 
neglected, typically because external sources are fixed). 

(iii) Quasisteady evolutions when both inertia and dissipation have to be included 
(system diffuses towards the critical state or towards a new steady state if dictated 
by the external sources). 

The rest of this paper will be focused on several examples in which the properties 
of type-II superconductors close to the critical state are calculated based on the theory 
issued in this section. The presentation will be organized according to the three cases 
described above. 

3. Numerical applications 

As it was explained in SecEj the core of our theoretical proposal is the existence 
of a dissipation function J-" that one may write down as a quadratic expression of the 
macroscopic current density vector components. Recall that both from the mathematical 
point of view, and also as concerns the physical background, J-" should be defined as 
a certain distance d. Actually, d is a measure of the separation between the operation 
point (given by the value of the current density J) and a certain equilibrium value Jc 
that depends on J itself. The specific form of the function J-" depends on the critical 
current region boundary Fc, as we have explicitly shown in Fig|2] for the case of an 
elliptic behavior. 

In this section, we present a number of examples in which the above ideas have 
been applied to various practical configurations. Let us start by writing down some 
useful equations. Just for clarity, our statements refer to the elliptical region in Figl2l 
but generalization is apparent. 

• Notice that one may classify a given value of J (in components (Jy, J±)) as lying 
inside F^ or not, by calculating the sign of the quantity 



FigED. 




(20) 
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i.e.: 5r(J) is either negative or positive when J is either within or beyond the 
contour of Fc (and null for the critical values Jc G r^). 

• The distance function for a given point {J\\,J±), also involving its image Jc = 
(J|*, J^) is given by 

d^=(j^^-j*y+(j^-Jl)\ (21) 

Here, J|* and have to be determined for the actual critical current law (region) 
under consideration. In the case of an elliptic model, the criterion of minimum 
distance between (Jy, J±) and the boundary leads to solving a quartic equation. 

• For a given value (Jy, J±), the dissipation function J^{3) may be written as 

^=iper(J)rf2(J) (22) 

with p the material resistivity, and d the distance function referred above. Gr 
stands for a step function whose value is zero within and one over the outside. 
A useful representation is 

^ , , 1 + tanh[Mr(J)l , . ^ 

0r(J) = ^^^^,A;»1. (23) 

Recall that for the rather usual experimental configuration in which J = (the 
components of J parallel to the magnetic field are zero), the above formulation may 
be cast as follows (mind the thick lines in Figl2]): 

Sr = {J±/JcL? - 1 



J? 



Jc_L if J_L > Jc± 

— Jc± if J± < —Jc± 



.2 _ ) (-'^l - Jc±) if J± > Jc± .X 

+ Jc±r if J± < -Jc± ^ ' 

These relations mean that, for a vast number of experimental setups, J-" may be built 
by composing two parabolae and a step function as illustrated in FigJSl The examples 
supplied below are based on this assumption. For the readers' sake, we eventually 
provide a useful discrete form of Eq. lfTSll that is based upon the transformation of 
the electromagnetic problem in terms of potentials [16] and on identifying the set of 
elementary circuits related to the problem's symmetry, viz. 

^ / ||Hn+i-Hn||W / AtJ^dV (25) 

2 Jr3 



Vol 



2 



Here means the vector potential component along the current lines, related to the 
applied magnetic field, Mij denotes the inductance between elementary currents Jj, Ij 
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Figure 5. (Color online) Detail of a one-dimensional dissipation function for type-II 
superconduetors. J- is built upon multiplying two parabolae by the step functions 
depicted below and adding them. Note that the horizontal axis, i.e.: J is common to 
boh plots. 



flowing along the speciflc circuits of the problem, and the tilded quantities concern 
previous time layer. At is the incremental time of our calculation, Ri stands for the 
incremental resistance of the i-th circuit and /c, its critical current. Application to 
different problems will imply different expressions for the matrix elements Mij, the 
applied vector potential, the resistance and the critical currents. 

Just for completeness, we mention that other possibilities for the dissipation 
function may be explored as it will be shown below. In fact, the customary power- 
law expression for the E{J) law will be also investigated in this paper, through the 
related dissipation function 

f AtJ^dV ^FoAty^lr-^ (26) 

ivol . \hiJ 

3.1. Metastatic equilibrium: the critical states 

To start with, we present an application of minimizing the expression in Eq. (l25l) in 
a genuine "critical" situation, i.e.: the dissipation term is not included exphcitly, but 
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Figure 6. (Color online) (Sheet) current density profiles in the equilibrium (critical) 
states that are induced in a superconducting disk by an increasing magnetic field 
along its axis (given by 2Ho/Jcd — 0.5,1,1.5,2,2.5,3,3.5 from the outermost to the 
innermost curve). Lines correspond to the analytical solution in Ref. [17j (see text) 
and symbols to our numerical results. 



replaced by the constraint J < J^. in the minimization of the first term. Physically, this 
corresponds to the conventional critical state problems in which overcritical excursions 
are neglected. We have chosen a case of interest, that has been solved analytically, thus 
allowing a straightforward comparison to our numerical results. In connection with the 
widely used configuration of a fiat sample in a perpendicular magnetic field, we apply 
our method to the circular disk geometry. According to Mikheenko and Kuzovlev |17j . 
the sheet current density distribution in a superconducting disk subject to a uniform 
perpendicular field of intensity Hq is given by 



with a = R/ cosh {2Hq/ Jed), d stands for the (small) thickness of the disk and R for its 
radius. Here, J* denotes the so-called sheet current that averages the current density 
over the thickness, i.e.: J*(r) = J J{r,z)dz. 

Figj6] displays the comparison of our results and those obtained from the above 
equation. We stress the fact that numerics produce the correct solution, even for 
situations in which the variable does not reach the true critical value. Recall that, 
physically, the strict inequality (J*(r) < J*) relates to the meaning of the sheet current 
as an average over the sample thickness. This situation occurs when a non-penetrated 
core still exists somewhere within the sample, i.e.: at some points {rQ,z) one has 
J(ro, z) = Jc for the peripheral values of z and J(ro, z) = in the central core. 




r < a 




(27) 
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3.2. Steady states: the voltage criterion 



Below, we analyze the opposite limit of the general conduction problem with type- 
II superconductors. The steady situation in which a dissipative state is maintained by 
the external action (current source) will be analyzed. Thus, based on the considerations 
introduced in Sec J2.2.3l for the law E = VjJ-", we concentrate on recently reported 
voltage-current experiments |5], specifically designed to investigate the complex E(J) 
law and underlying critical boundary Fc in High-Tc superconductors. Owing to the 
smart design of the experiment, a direct analytical study is allowed. 

First, we recall that the results in Ref. [5] show that the elliptic model for Fc gives 
an excellent fit to measurements. The authors have also obtained the polar angular 
dependence of the angle between E and J , say (3 (9) (see Fig|3] for definitions and, in 
our case, consider Eq. (fT6|) in order to obtain the direction of E relative to AJ). In 
their experimental configuration J is fixed along the film direction, while B is applied 
at different angles 6. The criterion for reaching and exceeding the critical state comes 
from a fixed value of E parallel to the direction of J {E\\ = E ■ J = Eq). 

From the elliptical curve of critical currents 



^Jc\\J \Jc±. 

we get the parametric representation 

'cos(a)\ ^ / sin(a) 



(28) 



Jria? 



Jc\\ J V Jc± 



(29) 



in terms of which the normal vector is given by (cos(q;)/ J^y, sin(a)/J^j_). Recalling the 
notation introduced in Sec j2.2.3[ fi = (cos(</)), sin(0)), we have 7^tan(0) = tan(a), with 
the anisotropy ratio 7 = Jc±/ Jc\\- 

Now, the ratio E±/E\\ for small dissipation (A J 0) is basically that of the electric 
field components at Fc, where E _L Fc {/3 = (p — 9), .i.e.: 



E±/E\\ = tan 



tan(^) - tan(0) tan{e){-f'^ 



' 1 + tan(0) tan(^) 72 + tan(^)2 

= ^c{i,e) (30) 

because in this limit a = 9. 

In order to quantify the influence of the dissipation parameters in the electric field 
ratio, we proceed by considering a first order correction in terms of A J. In the simplest 
case of constant p (see Sec J2.2.3]) this maintains E || n, but now a ^ 9. From the 
trigonometric relations in the triangle of sides J, Jc and A J and the condition E\\ = Eq 
we get ^ — a ^ ^cEo/{pJc), and then 
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74 + tan2(^) 

= $Aj(7,^,5) (31) 

where 6 = Eq/{pJc±) is the small parameter in the expansion. 

Eventually, we consider the most general case in which anisotropic resistivity is 
allowed. Let us assume the form p = cos^(^) + p_Lsin^(6'), and = P\\/p± the 
anisotropy ratio. Starting with Eq. (fT6l) and after some algebra we get 



(r2 + tan2(e))(74 + tan2(^))i/2 
= ^,r{l,0,6,r) (32) 

with $^ = (1 - r^) tan(6')/(r2 + tan2(6')). 

Notice that, considering that the anisotropy ratio 7 is given, our theory includes 
two parameters, 6 and r^, that may be used to fit the experiment. Ultimately this will 
either allow to obtain the resistivity coefficients or (if they are known by other means) 
check for consistency. 

Our theory has been used to analyze the experimental information in Ref . [5] • FigJT] 
shows our predictions as compared to the available data. Basically, we plot the angular 
dependence of the ratio E±/E\\ as given by the different levels of approximation in 
Eqs. fl30|31ll32l) . i.e.: $c; ^aj and $cr- Our theoretical curves are obtained by using 
7 = 0.34 the value that is reported by the authors as a direct information from their 
experiments. Several facts are of mention: 

• Experimental data have been "recalibrated", multiplying by a factor of 0.714 

• If one takes for granted that the value of 7 is correct, $c does not produce the best 
fit even when multiplying by constant factors (in the plot, we have the direct data 
and also multiplied by a = 0.714). On the contrary, allowing a constant factor 
between theory and experiment, the best results are given by $aj with 6 = 0.15. 

• the correction introduced by anisotropic resistivity ($«■) does not seem relevant to 
this experiment. Lacking direct information, we just see that changing its value 
from 3 to 1/3 only minor effects are observed. 

• If one analyzes the first derivative of $ and compares to the results obtained from 
the experimental data, some doubts arise, because does not seem to show the 
real shape, whereas <I>^ looks more realistic. 

From these issues, we state that the comparison between theory and experiment 
is reasonable. First, from the value 5 = 0.15 one obtains a resistivity: p ~ 10~^f2m, 
which is nor far from the typical values for flux flow resistivities at the experimental 
range (around 10~^Qm for YBCO films at T = 85/'^). Second, the calibration factor 
for the experimental data (C = 0, 714) is not unreasonable because the electric fields 
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Figure 7. (Color online) Plots of the electric field ratio Ey/Ez in the magnetotransprt 
experiment sketched in the inset. A transport current is directed along the z axis of 
a superconducting wire, that is also subjected to an applied magnetic field (induction 
at an angle 6). Electric fields are calculated for J ^ Jc according to Eas. (|30l3f I32p 
and plotted as lines. Symbols come from the experiment in Ref. [5]: YBCO film at 
T = 85K. Top left: comparison between our theoretical curves $c, ^aj, and ^cr for 
a couple of values of the anisotropy parameter (r = 0.3 & r = 3). Top right: $c is 
compared to the experimental data, either multiplied or not by a factor a. Bottom 
left: $AJ is compared to the experimental data. Bottom right: comparison of the 
derivatives of the function [Ey/ E2]{0) obtained from theory and experiment. See the 
text for details about parameter values in our simulations and scaling of data. 

were just obtained from the averaged expression voltage/ distance. According to the 
authors of Ref. [5] this procedure should be "refined" when (as it is the case in parallel 
configurations) the electric field structure along the sample is not homogeneous, a fact 
that is under investigation. 

3.3. Transient behavior: relaxation towards the critical state 

Finally, based on the minimization of the expression in the Eq. (l25l) . we put 
forward several cases that explicitly show the diffusion of electromagnetic fields in a 




Figure 8. (Color online) Illustration of the diffusion processes simulated in this work. 
Here, a step is applied to the transport current / along the X-axis of a thin strip 
(as shown in the insets). The circulating sheet current J(y) increases over the initial 
critical state profile (lower dashes) and then "relaxes" towards the subsequent critical 
state (upper dashes). Continuous lines correspond to the intermediate current density 
profiles, and arrows are used to indicate the time evolution of J at different parts of the 
sample. These relaxation curves are obtained subsequent to the step in the transport 
current (/—>/ + A/) and occur at equidistant time intervals of value 0.4 tq. 



hard superconductor. We will concentrate on the tape geometry (FigjH]) and analyze 
relaxation towards the critical state, when either transport current or magnetic field 
steps are applied. By considering different values of the bias characteristic period Tq as 
compared to the magnetic diffusion time constant Tp = fioL^/p (with L the sample's 
dimension), we will display different conditions in which critical states are (or aren't) 
realized after relaxation. It will be shown that the conditions studied may be applied 
to establish validity of the CSM approximation in AC experiments. 

The general features of our simulations are shown in FigJHl Notice that we obtain 
the current density profiles that occur between two equilibrium states when a step of 
transport current is introduced. Mainly, we have studied the infiuence of the ratio Tq/tp 
on the diffusion process within the quadratic dissipation function framework developed 
in the previous sections. This relates to the piecewise linear behavior of the E{J) law. 
Furthermore, for completeness, we have also analyzed the results for power-law relations 
as indicated above. 
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Figure 9. (Color online) Normalized sheet current (top) and magnetic field (bottom) 
profiles in a superconducting strip carrying an impressed transport current /, that 
is cycled as labelled in the curves (see text for the definition of units). Dashed 
lines correspond to the analytical solutions in Ref. |18j . and continuous lines to our 
numerical diffusion calculations. In this case, and for each current interval, we just 
plot the relaxation profiles obtained for a quadratic dissipation function at the time 
t ~ tq ~ 2.0 Tp, previous to the subsequent step in transport current. The panels to 
the left stand for the increasing branch of the applied current, while the right panels 
display the decreasing branch profiles. 



3.3.1. Superconducting strips with transport current 

A number of situations related to the experimental setup sketched in Fig|8] have 
been studied, but only a representative set of results are shown below (FigsJQl and ITOl). 
To start with, the choice of the long strip (tape) geometry allows a straightforward 
comparison to analytical results for the related limiting critical states. Thus, the thick 
dashed lines in that figure (critical states) have been obtained from the well-known 
expressions for the above defined sheet current [18] 



2 J'' la^ — W' 

h<\y\<a, 

with a the tape "half-thickness" and h = a^J{\ — P / 1"^^.^, /max = 2aJ^ = 
27raHc. In our case, normalized units have been used for all the physical quantities: 
J/ Jc, H/Hc, I / Iraa.x,y /d- The magnetic field profiles have been obtained by integration 
of the current density. Analytical expressions derived from Eg. (1331) may be found in 
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Figure 10. (Color online) Same as FigEl but now for tq/tp = 0.2. 

Ref. [IH]. As for the results of this paper, straightforward numerical integration of the 
current density data was performed. 

FigJH shows the penetration profiles for a transport current that increases from the 
virgin state until the value ///max = 0.98 is reached. Then, a negative ramp towards the 
value ///max = —0.98 was applied. The associated field penetration profiles are shown 
in the lower panels. Just for the readers' sake, we recall that expressions for critical 
state curves in the negative ramp of transport current may be obtained from Eq. (|33|) by 
applying linear superposition [18]. For the numerical curves in this plot, the calculated 
relaxation takes place under the condition tq/tp = 2 and here we just show the profiles 
previous to the application of the following step (evolution takes place in the fashion 
shown in FigJH]). Notice that for the increasing branch, a high degree of coincidence 
between the "relaxed" profiles and the exact critical states displays. It is only for the 
nearly penetrated sample that relaxation is noticeably incomplete at t = tq = 2.0 Tp. A 
mismatch appears in the central part of the tape, that is inherited by the decreasing 
branch of the cycle, i.e.: when new steps are applied before reaching the equilibrium 
profile. As a counterpart, a reasonable coincidence is always observed at the sample 
edges. Apparently, this behavior relates to the concept that flux penetrates from the 
surface, and that dissipation is an integral over the sample, thus giving place to a most 
effective relaxation when a smaller region is involved. In other words, relaxation leads 
to equilibrium quicker close to the surface than in the bulk. 

FigJTO] has been obtained in a similar fashion, but under the condition tq/tp = 0.2. 
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As expected, clearly incomplete relaxation is obtained. In practice, this would lead to 
magnetization profiles with absolute values beyond the critical state limit, and also to 
higher AC losses, corresponding to the noticeable excursions of the transport current 
density towards the region J > J^. 

We conclude this part by mentioning that results obtained for the power law relation 
are rather similar to the previously described behavior, when the respective choices 
m = 100 and m = 10 are done (to be compared to the cases Tq/tp = 2 and To/vp = 0.2 
respectively). 

3.3.2. Superconducting strips with applied magnetic field 

The diffusing current and field profiles in a thin strip in a perpendicular field Hq 
and with zero transport current have been calculated under a wide set of conditions. 
Figs {11] and [12] present the main features observed in our simulations. In this case, the 




Figure 11. (Color online) Normalized sheet current (top) and magnetic field 
(bottom) profiles for a superconducting strip in an impressed perpendicular magnetic 
field Hq. To the left, the depicted profiles correspond to the values Hq/Hc = 
0.16,0.32,0.48,0.64,0.8. To the right, the applied field is cycled, i.e.: Hq/Hc = 
0.8, 0.48, 0.16, —0.16, —0.48, —0.8. Dashed lines correspond to the analytical solutions 
in Ref. |18| . and continuous lines to our numerical diffusion calculations. In this case, 
we plot the profiles obtained for a quadratic dissipation function with tq/tp = 0.2 just 
previous to the subsequent step in transport current. Notice that the upper left plot 
also includes the results for tq/tp = 2.0, that are omitted in the rest to avoid confusion, 
and virtually coincide with the analytical lines. Labeling is bypassed in some apparent 
cases. 
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seed for the analytical evaluations is the expression 



2J' 



tan 



cy 



TT 



\y\ < b 

b < \y\ < a , 



(34) 



where c = tanh(ifo/-f^c)- Magnetic field profiles may be obtained by integration and 
negative ramp equations by linear superposition [18]. Also as before, dashed lines stand 
for the analytical results of equilibrium profiles, and continuous curves for our numerical 
calculations of finite time (incomplete) relaxation processes. We only display the profiles 
corresponding to the diffusion step just previous to the change of external condition. 

Fig{TT] shows the results for the quadratic dissipation function, and Fig{T2] for 
the power-law relation. It is noticeable that, in both cases, a remarkable degree of 
coincidence between the relaxed profiles and the eventual critical state solution occurs 
for Tq/tp > 2. Nevertheless, important differences appear for the case Tq/tp ^ 0.2. As 
related to such differences, an important feature has been observed. We call the readers' 
attention that in FigjlT] J*(t/) displays a basically linear penetration profile close to the 
sample's edge, whereas a kind of square root trend is observed in FigJT2l This may be 
explained as follows. According to Faraday's law, one has dyE^ = dtB^. Then, for small 
variations dtBz will be practically constant in space and so that E^ ~ (dtHo) y, which 
implies a linear behavior of J^{y) in the case of a piece- wise linear E{J) law, and an 
inverse power law when E and then, ~ y~^. 




Figure 12. (Color online) Same as FigfTTl but for the power law model with m = 10 
and m = 100. 
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Figure 13. (Color online) Normalized sheet current (top) and magnetic field 
(bottom) profiles for a "field-like" state in superconducting strip with both transport 
current and perpendicular magnetic field applied with constant ratio. In normalized 
units, the different steps correspond to Hq/Hc = 0.08, 0.16, ...0.72 and I/Ima^ = 
0.025, 0.05, ...0.225. Dashed lines correspond to the analytical solutions in Ref. jl8j . 
and continuous lines to our numerical diffusion calculations. In this case, we plot the 
profiles obtained for a quadratic dissipation function with tq/tp = 20. 



3.3.3. Superconducting strips with transport current and applied field 

Just for completeness, we will also provide an example of relaxation towards the 
critical state when both a transport current and a magnetic field are applied to the 
superconducting strip. Analytical expressions for the case of synchronous ramps of 
current and field have also been provided in Ref. [H] . Thus, the solution is built through 
the "generating function" 

b < y < a 

j{y,a,b)={ ^cot-'"—^ \y\<b (35) 

— oo < y < —b , 




where p = (y^ — b'^){a'^ — b"^). Explicitly, one has 

.r{y) = J' [j{y + w,a + w,b) + pj{-y -w,a-w, b)] (36) 
with the definition w = (///max)/ tanh(//o///c). 
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FigJTSl shows the comparison of our numerical results and the analytical calculation 
described above. Following the nomenclature in Ref. [18] we have concentrated on a 
so-called field-like situation given by r = I/27iaHo= 0.32. 

In this case, we have considered a high relaxation ratio Tq/tp = 20, and a practical 
coincidence with the critical state profiles is observed. 

4. Conclusions 

The motivation of this work has been to explore the possibility of extending 
the critical state concept so as to include the effects of flux flow resistance in a 
general sense, i.e.: allowing the possibility of overcritical current density in several 
dimensions. More specifically, we have studied the case of homogeneous type-II 
samples, in which anisotropy is introduced by the direction of the local magnetic 
field. However, generalization to other cases, such as microstructure induced anisotropy 
are straightforwardly dealt with. Conceptually, we have chosen a macroscopic 
(thermodynamic) point of view, that avoids the explicit consideration of the underlying 
vortex physics. In brief, our theory relies on two facts: (i) the experimental evidence that 
in type-II superconductors dissipativeless currents are allowed when the components 
of J either parallel or perpendicular to some specific direction do not exceed specific 
thresholds, and (ii) the resistive transition has to occur according to the laws of entropy 
production. These considerations are formulated in a geometrical language through the 
definition of the so-called dissipation function, that takes values over the space of current 
densities, i.e.: J^{J). T goes to zero within the so-called critical region A (generalization 
of the one dimensional condition J < Jc) and is a positive definite quadratic form 
beyond. The main role of this function relates to the issue of the material law E(J) 
that appears by just imposing consistency and uniqueness in the physical quantities. In 
fact, one has E = Vj-F. 

The relevant approximations of the theory may be stated in terms of time constants. 
As related to the typical time scale of the macroscopic experiments (i.e.: IHz < cu < 
IkHz), the following scale applies: (i) charge recombination processes are assumed to 
occur instantly (magnetoquasistatic approximation), (ii) magnetic diffusion occurs very 
fast in the Critical State limit and may be modulated by a time constant Tp oc 1/p in 
our theory, and finally (iii) thermal activation and relaxation to the true equilibrium 
(J = 0) occurs very slowly and may be neglected. 

From the practical point of view, and taking advantage of the variational 
interpretation of the electromagnetic problem [19], we put forward a minimization 
statement that gives way to the numerical form of our theory. Eq.f l25p displays the 
function that is minimized, and is the central result of this article. One can identify 
two basic contributions that relate to the physics of the problem: (i) the inertial terms 
that account for the reversible energy storage, and (ii) the energy dissipation term that 
includes T . Inspired by this interpretation, we have worked out a number of examples 



Electromagnetics close beyond the critical state: thermodynamic prospect 



30 



that illustrate the application of the theory either straightforwardly to the standard 
critical state problem, to steady states in which permanent dissipation is forced by some 
external action, or to the transient (diffusive) processes that occur in between successive 
critical states, or between steady states as dictated by modified external sources. 

Based on the comparison of our results to other theoretical works on the critical 
state, we conclude that the complex non-linear diffusion processes that take place 
according to the theory, converge to the previously reported critical states solutions 
for a given set of external conditions. Excellent agreement is observed when the system 
is allowed to relax within the typical excitation period. 

Related to the analysis of dissipative steady states, we stress the fact that our theory 
unifies the standard CSM framework and the results of current- volt age techniques, used 
to derive the critical current parameters. In this sense, we give a number of relations 
(Eqs.(l30l) to ( 132|) ) that allow to analyze the multicomponent E(J) measurements 
designed for characterizing the critical current behavior of High-Tc superconductors [5] • 
Along this line, we have compared our predictions to those experiments and can conclude 
that, although a reasonable agreement is obtained, more experimental information would 
be desirable in order to choose between one theory or another for the overcritical states. 

Although we state that a quadratic dissipation function (and thus a piecewise linear 
E( J) law) is the more judicious choice for investigating the behavior close to the critical 
state, the rather extended use of a power-law relation has also been checked against our 
results with good degree of coincidence. Only, second order qualitative differences can 
be observed. Notice, in passing, that big discrepancies with the critical state are out 
of reach of our theory, but also of experiments because holding current densities much 
above the critical value would result in a blow-up of the sample. 

Further work along the lines of this paper entails the application of the theory 
to obtain the relaxation profiles in higher dimensional systems, such as structurally 
anisotropic materials and problems with multicomponent magnetic fields with non 
homogeneous electric field profiles. 
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